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ABSTRACT 

We present the data processing and analysis techniques we are using to determine structural and 
photometric properties of galaxies in our Gemini/HST Galaxy Cluster Project sample. The goal of this 
study is to understand cluster galaxy evolution in terms of scaling relations and structural properties 
of cluster galaxies at redshifts 0.15 < z < 1.0. To derive parameters such as total magnitude, half-light 
radius, effective surface brightness, and Sersic n, we fit r 1 / 4 law and Sersic function 2-D surface brightness 
profiles to each of the galaxies in our sample. Using simulated galaxies, we test how the assumed profile 
affects the derived parameters and how the uncertainties affect our Fundamental Plane results. We find 
that while fitting galaxies which have Sersic index n < 4 with r 1 / 4 law profiles systematically overes- 
timates the galaxy radius and flux, the combination of profile parameters that enter the Fundamental 
Plane has uncertainties that are small. Average systematic offsets and associated random uncertainties 
in magnitude and logr e for n > 2 galaxies fitted with r 1 / 4 law profiles are —0.1 ± 0.3 and 0.1 ± 0.2 
respectively. The combination of effective radius and surface brightness, logr e — /31og(/) e , that enters 
the Fundamental Plane produces offsets smaller than —0.02 ±0.10. This systematic error is insignificant 
and independent of galaxy magnitude or size. A catalog of photometry and surface brightness profile 
parameters is presented for three of the clusters in our sample, RX J0142. 0+2131, RX J0152. 7-1357, and 
RX J1226.9+3332 at redshifts 0.28, 0.83, and 0.89 respectively. 

Subject headings: galaxy clusters: individual RX J0142. 0+2131, RX J0152. 7-1357, RX J1226.9+3332 - 
methods: data analysis - galaxies: fundamental parameters 
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1. INTRODUCTION 

Theoretical hierarchical models of galaxy formation pre- 
dict that numerous small halos, primordial dwarf galaxies, 
collapsed early on in the history of the universe. These 
merged to form ever larger halos with deeper potentials 
and greater baryonic mass accumulation. While these 
ACDM models are successful at explaining structure for- 
mation on large scales, the agree ment with observation s 
is worse on the level of galaxies ([Conselice et al.1 120071) . 
Understanding the formation of primordial galaxies and 
how galaxies came to have the stellar populations, masses, 
and structural properties observed in the local universe is 
fundamental to our understanding of cosmology. 

Generally, the physical processes a galaxy undergoes 
during its evolution depend primarily on only two fac- 
tors: mass and environment. While mass may affect 
galaxy structural and dynamical properties and star for- 
mation activity in predictable ways, the local environment 
of the galaxy is the source of complex evolutionary pro- 
cesses such as galaxy-galaxy merging, harassment, and 
stripping. Despite this, tight scaling relations of galaxy 
properties exist and must be reconciled with potentially 
different evolutionary paths. Such scaling relations in- 



clude the Tully-Fisher relation bet ween disk galaxy tota l 
magnitude and rotational velocity (jTullv fc Fisher!lT977t ). 
the three-parameter Fundamental Plane (FP) relationship 
between early type galaxy velocity dispersi on, effective ra- 
dius, and effective surface brig htness (e.g. iDressler et alJ 
ll987tlDTorgovski fc Davislll987f ). the red sequence for early 
type galaxies, and relationships between absorption line 
strengths and central velocity dispersion. How these tight 
relations evolve with redshift or vary with galaxy mass 
places strong constraints on galaxy formation and evolu- 
tion models. 

Observational studies of cluster galaxies have found 
that the most massive galaxies are composed mainly of 
old stellar populations suggesting that massive early type 



|j0rgensen et al. 


1999: 


l2003t iMei et al. 


2006) 



[Traeer et al.l 120001: iBlakeslee et all 
. There is little evidence for recent 
star formation in these quiescent galaxies. They form a red 
sequ ence with little scatter in a col or-magnitude diagram 
(e.g. ISandage fc Visvanathanl |T978) and they obey tight 
scaling rela tions in their kinema t ic and structural prop - 
erties (e.g. IDressler etaLl Il987t Ugrgensen et al.l 1 19961 ). 
These all suggest an early epoch of star formation in a 
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homogeneously old galaxy population. Although mas- 
sive galaxies exhibit pure passive evolution since z ~ 2 
evolving only in luminosity and color as stellar popula- 
tions age, observations have also shown that clusters and 
lower mass cluster galaxies continue to evolve at interme- 
diate redshifts. The fraction of blue, star forming galax- 
ies increases with re dshift (e.g. iButcher fc Oemleii Il984t 
Ellingson ct al. 2001) while larger fractions of spiral to SO 
gala xies are observed in clusters at intermediate redshifts 
(e.g. IDressler et aT]|1997l) . 

Meanwhile hierarchical structure formation models pre- 
dict that galaxies form from the collective mergers of 
smaller galaxies. These mergers are expected to continue 
with high frequen cy through intermediate redshifts (e.g. 
iBaugh et al.lll996h . Semi- analytical models predict a char- 
ac teristic mass a factor of 3 lower than found at z = 
1 (|Poggiantil 120041 ) with simulations underpredicting the 
numbers and mass densi ties of the most mass ive galaxies 
by a factor of over 100 (jConselice et al.l [20071 ). all imply- 
ing that massive galaxies undergo a faster build-up than 
predicted. Dynamical models also often show a greater 
range in e arly type galaxy structu re and kinematics than 
observed (|de Zeeuw fc Franxlll991l ). 

In order to reconcile apparently contradictory observa- 
tions with formation and evolution models, we must study 
galaxy populations over a wide range of redshifts and for a 
wide range in galaxy mass. Previous studies have tended 
to investigate larg e samples of galaxies within a single c lus- 
ter or epoch (e.g. iFritz et all [20051 ; ISmith et al.ll2004D. or 
fewer galaxies over a range of redshifts (e.g. iTreu et ail 
2002]), and often only the brightest, most massive galax- 
ies. Through the Gemini/HST Galaxy Cluster Project, 
we seek to better understand the processes driving cluster 
galaxy evolution by studying scaling relations as a func- 
tion of mass and environment since z = 1.0, at a time when 
the universe was only half its current age. We have there- 
fore obtained data for a large sample of galaxies, which 
includes a sufficient number of galaxies at each redshift to 
accurately measure scaling relations and a sufficient num- 
ber of redshifts to measu re trends in the scaling rela tions 
as a function of redshift (jJprgensen et al.|[2006l l2007l ). At 
each redshift our sample spans a wide range in galaxy lu- 
minosity in order to investigate the role of galaxy mass. To 
discriminate between evolutionary effects due to environ- 
ment and due to galaxy mass we have obtained consistent 
radial coverage for each cluster. 

In a previous paper we hay e described the project and 
sample ( Jpr gensen et al.ll2005l) . Our sample of galaxy clus- 
ters consists of 15 X-ray selected massive clusters with red- 
shifts ranging from 0.15 < z < 1. For each of the clusters 
in our sample, we have obtained Gemini/GMOS spectra 
for ~ 30 cluster members with a wide range of luminosi- 
ties. The targetted objects are chosen independently of 
morphology since there is evidence that morphologies may 
continue to evolve even to t he present epoch and w e wish 
to avoid "progenitor bias" (|van Dokkum fc Franxll200"ll) . 
From the spectra, we obtain redshifts, velocity disper- 
sions, and line index measurements. Using HST imag- 
ing data, we measure structural parameters and surface 
brightness profiles. With this large sample, we are study- 
ing scaling relations such as line index strengths and the 
FP as a function of redshift. Line index measurements 



probe chemical enrichment providing insight into the star 
formation histories of the galaxies. The FP is examined 
in terms of mass and mass-to-light (M/L) ratios in order 
to probe both assembly histories and luminosity evolu- 
tion. We are using structural parameters for quantitative 
measurements of morphology to study morphological evo- 
lution. We have published t he results on line inde x re- 
lations for RX J01 52. 7-1357 fj0rgensen et al.l l2005h and 
RX J0142.0+2 131 dBarr et al.f[2005f) . and the FP for RX 
J0142.0+2131 dBarr et al.ll2006fl and RX J0152. 7-1357 and 
RX J1226. 9+3332 (|Jc3rgensen et al.ll2006l . [200l . 

The FP relates surface brightness, effective radius, 
and velo city dispersion in a tight and well established 
relation (IDressler et"all Il987t iDiorgovski fc Davisl Il987t 
Uprgensen et al.lll996f) 

log r e = a log a + /31og(/) e +7. (1) 
It can also be described as a relation between M/L and a 
or as 

log(M/L) =£logM + 7' (2) 

which makes the FP a powerful tool to study the star 
formation and assembly history in early type galaxies. 
A comparison of the FP for our high redshift sample to 
that of low redshift Coma cluster galaxies has revealed 
a number of interesting results. The two z ~ 0.85 clus- 
ters exhibit a steeper slope than the low redshift FP. We 
find this to be evidence for downsizing in which the lower 
mass galaxies have undergone more recent star forma- 
tion and are overluminous compared to their low-z coun- 
terparts (|j0rgensen et al.ll2006l l2007f) . The FP for RX 
J0142. 0+2131 at z=0.28, while having the same slope 
as the low z sample, displays a greater scatter, possible 
evidence for the galaxies having undergone rapid bursts 
of star formation during a cluster merger at z > 0.85 
(|Barr et al.ll2006h . 

Future papers will address the FP, line index scaling re- 
lations, and galaxy quantitative morphologies for all clus- 
ters in our sample in order to compare the star forma- 
tion and assembly histories over a range of redshifts. Here 
we present the data reduction and analysis techniques we 
have used to measure structural and photometric param- 
eters for three clusters in our sample, RX J0152. 7-1357, 
RX J1226. 9+3332, and RX J0142.0+2131. 

RX J0152. 7-1357 is a massive cluster at z=0.83 orig- 
inally discovered from ROSAT data. XMM-Newton and 
Chandra observations later showed the cluster to consist of 



two subclumps 


in the early stages of merging ( 


Jones et alj 


120041: iMaughan 


et al.H2003HGirardi et al.ll2005 


) . Much re- 



search has been done on th is cluster, including r ecent stud- 
ies of star formation rates (|Homeier et a l. 2005), morphol- 
ogy, and the color-magnitude diagram (jBlakeslee et al.l 
2006). RX J1226.9+3332 is a massive cluster at z=0.89 
and also the most X-ray luminous cluster known at such 
high redshift. It was discovered in the Wide Angle ROSAT 
Pointed Survey (|Ebeling et al.1 12001D and exhibits a re- 
laxed morphology. However, in a deep XXM-Newton and 
Chandra study of the X-ray mass analysis of this cluster, 
IMaughan et alj (|2007f ) find evidence for a recent or on- 
going merger event. RX J0142.0+2131, at z=0.28, was 
first identified as a massive cluster in bo th the Northern 
ROS AT All-Sky Galaxy Cluster Survey (|Bohringer et al.l 
120001) and the ROSAT extended Brightest Cluster Sam- 
ple (|Ebeling et aLll2000l ). Although it is a relatively poor 
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cluster, it displays a large clust er velocity dispe rsion, yet 
shows no signs of substructure (jBarr et al.ll2005l ). 

In this paper we describe our measurements of the phys- 
ical properties of the galaxies in our sample and analyze 
how uncertainties affect the derivation of the scaling rela- 
tions and overall results of our study. In Sections 2 and 
3 we describe our observations and our imaging data re- 
duction process. Catalogs of photometric and structural 
parameters for RX J0152.7-1357, RX J1226.9+3332, and 
RX J0142. 0+2131 are presented in Section 4. We discuss 
internal and external consistency of our results and impli- 
cations on the measurements of the FP in Section 5. A 
summary is presented in Section 6. We assume Hq = 70 
km s- 1 Mpc" 1 , TO = 0.3, and fi A = 0.7. 

2. OBSERVATIONS 

We have acquired ground-based imaging for the clus- 
ters in our sample using GMOS on the Gemini-North 
8m telescope. GMOS spectroscopy was obtained for 
40 — 50 galaxies in each cluster field, 20 — 30 of which 
turned out t o be cluster members. GMOS is described in 
iHook et al. ( 2004). Details of how the spectroscopic sam- 
ple was chosen and the observations and reduction of these 
data are described elsewhere (J0rgens en et al. in prep; 
iJprgensen et al.l 120051 : iBarr et all [2005) . For each of our 
clusters we have obtained new ACS imaging or use archival 
HST data to study the galaxy structural parameters with 
high resolution data having FWHM = 0".l. Data for the 
three clusters described in this paper were obtained using 
the ACS Wide Field Channel. The ACS/WFC detector 
consists of two 2048x4096 chips separated by a gap of 2.5 
arcsec with a mean pixel scale 0" .049 and a 3'. 36 x 3'. 36 
field of view (jSirianni et al.ll2005f h 

For RX J0152. 7-1357 (z=0.83), we used archive data 
from Program ID 9290 which includes 4 mosaicked ACS 
fields observed in 3 filters each, F850LP, F775W, and 
F625W, corresponding to the z', i', and r'— bands in our 
GMOS data. Total exposure times for each field were 
4800s in the i'—band and 4750s in the z' and r' — bands. 
These data contain all 29 of our spectroscopic sample clus- 
ter member galaxies. 

For RX J1226.9+3332 (z=0.89), ACS archival data 
from Program ID 9033 included four fields in two bands 
each, F606W and F814W, the latter equivalent to the 
Cousins /—band. Exposure times were 8 x 500s. These 
fields contain our original 25 sample galaxies used in 
Uorgensen et al.l (|2006l l2007h along with an additional 87 
galaxies from an extended sample. 

We obtained data for RX J0142.0+2131 (z=0.28) from 
Cycle 12 Program ID 9770 observed on UT 2003 November 
01 and 2004 July 03. Two fields were imaged in a single 
band, F775W. The total exposure time for each position 
was 4420s. 28 of 30 spectroscopic sample cluster member 
galaxies were covered by these fields. See Table [1] for a 
summary of all the ACS data used in this paper. 

3. DATA PROCESSING 



3.1. HST / ACS Initial Processing and photometry 

The three clusters discussed in this paper were used to 
establish the data reduction pipeline. Although the data 
are distributed already processed from STScI we chose to 
redo the reduction starting with the stacking of flattened 
images in order to achieve better cosmic ray removal. In 
some cases we also wished to mask out large reflections 
or saturated stars before stacking the images. We used 
MULTIDRIZ ZLE v. 2 3, an STScI distributed PyRAF 
taskQ (jKoekemoer et all |2002|) . which drizzles, distortion 
corrects, and stacks images all in one. 

To obtain better relative shifts between separate expo- 
sures at each pointing, we initially ran MULTIDRIZZLE 
to generate a median combined image for each pointing 
along with individually drizzled images shifted to the ref- 
erence frame of the first image in each set. The images at 
this point were aligned simply by using the header WCS. 
We then ran a subset of the Gemini IRAF packag^B task 
IMCOADD which first detects stars on the cosmic ray- 
free median image using DAOPHOT. It then locates these 
stars on the individual drizzled images and measures any 
additional relative shifts. This step created a shift file suit- 
able as input for MULTIDRIZZLE which was subsequently 
run a second time to create the final stacked, distortion 
corrected, cosmic ray cleaned image. A wrapper script 
was used to perform the necessary steps which included 
an option to allow manual sky subtraction, for use when 
bright stars or reflections effected the image. We chose not 
to perform any resampling during the drizzle and we used 
the default square drizzle kernel as no improvement was 
observed with other kernels. 

MULTIDRIZZLE output includes weight and context 
images. Each pixel within the context image is encoded 
with information about which specific individual images 
were used during the combining operation for that pixel 
(ie. were free from chip defects and CRs). We use 
the context image to generate a map which more gen- 
erally contains information about the number of pixels 
that were combined to create each stacked pixel value. 
This is used to calculate more rigorously the sigma noise 
map for each ACS image taking into ac count the corre- 
lated noise produced by the drizzling (jCasertano et al.l 
l2000f) . The noise at each pixel is calculated as a 2 = 
F A ((RN * VNcomb/gaineff) 2 + {N/gain e ff)) where RN 
is the readnoise, N com {, is the number of pixels used in the 
drizzling, gain e ^ f is the effective gain at each pixel incor- 
porating N com (,, and N denotes the counts in that pixel. 
F^ corrects for the correlate d noise and comes from ap- 
pendix 6 in iCasertano et al.l (|2000f ). It depends on the 
output pixel scale (in this case the fractional size com- 
pared to the original) and the drop size for the drizzling. 
In our case, scale and pixfrac are both 1.0 since we do not 
drizzle to smaller scales. F^ is also dependent on the size 
of the area considered since correlated noise scales from 
single pixel noise differently than uncorrelated noise. We 
did not take this area correction into account here. Thus, 
our noise map may incorrectly estimate the true photo- 
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metric errors. Although this would have an effect on the 
GALFIT error estimates on the individual parameters, we 
do not expect this to affect our results since we deter- 
mine parameter measurement uncertainties independently 
through galaxy simulations (Section [5]). 

The weight image uses weighting based on exposure 
times. To improve extraction of galaxies in later process- 
ing using SExtractor, we adjust this weight image using 
our n-pixel map described above to set a very low weight 
near image borders and in gap regions where few overlap- 
ping images produce greater noise. 

The zcropoint calibration in the AB system is pro- 
vided by STScI for the ACS camera in each filter used 
( Siri anni et al.|[2005h . We correct the WCS in the header 
of the stacked images, which were in most cases off by a 
shift of about 2" in RA and 1" in DEC, by using the IRAF 
tasks CCMAP and CCFIND along with coordinates for 
these galaxies previously obtain ed from our GMOS images 
relative to the USNO catalog (jMonet et al.lll998t ). Final 
coordinates are generally good to about 0.2" in both RA 
and DEC. 

3.2. SExtractor Object Detection and Photometry 

We used SExtractor (|Bertin fc Arnoutslll996f ) v.2.3.2 to 
detect all objects in the images having at least 9 contigu- 
ous pixels greater than 2a above the sky noise. The weight 
image described above was used to minimize the number 
of spurious noise detections predominently found near im- 
age borders and, in some cases, chip gap regions of the 
combined dithered images. 

SExtractor was run in association matching mode to 
specifically recover our spectroscopic sample galaxies. It 
was then run a second time to detect all other galaxies in 
the field, using dual image mode in cases where fields had 
been observed in multiple bands. In this mode, we use the 
band closest to the i'— band for the initial detection of the 
objects. After using the IRAF task imshift to precisely 
align images in other bands to the i 1 — band image, we run 
SExtractor to perform photometry in these other bands 
using the i'— band determined apertures. This produces 
matched catalogs of objects in all bands. 

After detecting all objects in a first pass through 
the data, SExtractor makes a second pass to deblcnd 
merged objects. We ran a number of tests to determine 
the best set of SExtractor deblending parameters, DE- 
BLEND_NTHRESH (sets the number of deblending lev- 
els) and DEBLEND.MINCONT (minimum fraction of pix- 
els a branch must have to be considered a separate ob- 
ject), which would correctly split and detect objects in 
our fields. Aperture check-images were used to determine 
whether objects had been detected correctly. For our spec- 
troscopic sample galaxies, the best set was found to be 
DEBLEND.NTHRESH = 32 and DEBLENDJVUNCONT 
= 0.01 as these parameters best recovered galaxies in their 
entirety without excess splitting while still correctly sep- 
arating neighboring objects. While this worked for most 
program galaxies in all clusters, there were cases of galaxies 
with a lot of structure that were split too much, or galaxies 
with foreground stars which required greater deblending. 
For these galaxies we adjusted the DEBLEND_MINCONT 
slightly lower or higher as necessary until that individual 
galaxy was verified by eye to have been detected correctly. 



For the non-sample galaxies in the images, many consist- 
ing of late type, irregular galaxies, we used a compromise 
set of parameters of DEBLEND_NTHRESH = 16 and DE- 
BLEND.MINCONT = 0.01. This minimized oversplitting 
while still separating most neighboring objects. 

The SExtractor output, which we use primarily as input 
for more rigorous galaxy profile fitting, included various 
magnitude (and flux) determinations, fractional pixel cen- 
troids, object class, ellipticity, position angle, and various 
moment and size measurements. We take the parameter 
nibest, a combination of adapted aperture and corrected 
isophotal magnitudes, as the estimate for total magnitude. 
Measured parameters for objects detected multiple times 
in overlapping image frames were averaged together. 

3.3. GALFIT Surface Photometry 

To perform surface brightness profile fitting, we used 
GALFIT , a two-dimens ional profile fitting program writ- 
ten bv lPeng et "all (|2002t l. GALFIT has the advantage over 
other similar profile fitting codes because it simultaneously 
fits neighboring galaxies. This improves the fit for r 1 / 4 law 
galaxies which contain a sig nificant amount of lig ht in the 
outer wings of the profile (|Haussler et al.l 120071 ). GAL- 
FIT provides the user with a choice of analytic surface 
brightness profile functi ons w ith which to fit the galaxies 
including Sersic (|Sersidll968h . de Vaucouleurs (r 1 / 4 law), 
exponential disk, bulge/disk deconvolution, and psf. 

The user must provide initial guesses for the profile pa- 
rameters and a PSF template used to deconvolve the image 
PSF during fitting. We use the SExtractor output m;, est 
for total magnitude, x, y, ellipticity, theta (corrected to 
GALFIT orientation) , and flux radius (set within SExtrac- 
tor to be an estimate of the half-light radius and converted 
to a e// for GALFIT) as input for GALFIT. Tests were 
made using simulated galaxies to determine the most ap- 
propriate PSF model and size. See Sect ion ISTTl for more de- 
tails of this testing process. The chosen P SFs were created 
by first running Tiny Tim (|KristJll995[ ) to generate raw 
distorted ACS model PSFs for a grid of locations spread 
uniformly over the two ACS chips. These were then added 
to a set of blank ACS images with the appropriate shifts 
and multidrizzled to produce PSFs combined in the same 
manner as the real data. These PSFs were re-extracted 
from the rectified frame as individual PSF templates for 
use with GALFIT. We use a 1-time sampled 9" PSF for 
the GALFIT fitting procedure. 

We fit all galaxies with both r 1 / 4 law and Sersic (|Sersid 
1968) profiles. The Sersic profile has the form 



E(r) = S e exp[- K ((r/r e ) 1 /"-l)] 



(3) 



where r e is the half-light radius, £ is the surface bright- 
ness, and n is the Sersic index, k is coupled to n to ensure 
that half of the total flux lies within r e . The total flux can 
be calculated by integrating out to r = oo: 



F tot = 2wr 2 e i: e e K n K - 2n T(2n)q/R(c) 



(4) 



R(c) = Ac + 2)/(4/3(l/(c + 2), 1 + l/(c + 2))) (5) 

(|Pend 120011 ) where q is the axial ratio, b/a, c is the 
diskiness-boxiness parameter, and (3 is the Beta function. 
The advantage of fitting with a Sersic profile is that the ex- 
act form of the profile need not be known a priori since the 
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index, n, can vary to fit the full range of possible galaxy 
profiles from n = f exponential disks to n = 4 ellipti- 
cals/spheroidal components. 

Since ACS images are large, we used a wrapper script 
which first extracted a smaller panel for each galaxy from 
the image, with a panel size determined to be ~ 25xr e of 
the galaxy to be fitted, having a minimum size of 250 pixels 
and maximum of 1000 pixels on a side. Because light from 
the extended profiles of neighboring objects will influence 
the profile fitting of the primary galaxy, all galaxies within 
the panel were either fitted or masked. Masking was nec- 
essary to limit the number of galaxies fit and speed up the 
execution of the fitting. We therefore simply masked out 
all objects detected by SExtractor more than four magni- 
tudes below the mean sample galaxy brightness. All other 
objects in the panel were fitted with Sersic profiles. These 
image sections were also checked by eye to ensure that no 
objects had been missed by SExtractor. 

The wrapper script fit all cluster sample galaxies with 
both a Sersic and a de Vaucouleurs (r 1 / 4 law) profile. For 
each profile we have a minimum of 6 free parameters, fit- 
ting for total magnitude, effective radius, ellipticity, po- 
sition angle, and x,y position. In the case of the Sersic 
profile, the Sersic n was set initially at 1.5 and allowed 
to vary. After fits with both Sersic and r 1 / 4 law profiles 
were performed, any object best fit with n > 3.0, was re- 
fit with a Sersic profile using n=4.0 along with the output 
from the r 1 / 4 fitting as initial guesses. All parameters were 
allowed to vary. We included this iteration because the ini- 
tial guesses tended to influence the fit and SExtractor did 
not always produce the most accurate estimates for input. 
Neighbors were fit simultaneously with Sersic profiles while 
the sky value was held constant at the value provided by 
SExtractor. Best fits are determined by minimizing the 
\ 2 residual, adjusting all free parameters simultaneously. 

Final parameters produced by GALFIT include total 
magnitude, effective radius, axial ratio, position angle, x 
and y centroids, Sersic n, and a xt value of the fit, along 
with associated uncertainties for all the derived parame- 
ters. Output also consisted of a multiple extension image 
including the original panel, a model image, and the resid- 
uals from the fit. We display several examples of r 1 / 4 law 
fits in Figures [T] - [3] Larger residuals are evident in galax- 
ies with best fit Sersic index n closer to 1 and in those 
exhibiting features beyond a smooth profile. Final xt va l~ 
ues typically indicated good fits. Median values of \1 f° r 
fits with both profiles in all three clusters were ~ 1.4 with 
maximum values for a few galaxies between 4 — 8. These 
were cases of galaxies with structures that were not well 
modeled with single profiles or of galaxies with problem- 
atic neighbors. 

GALFIT quoted uncertainties generally underestimate 
the true uncertainties of the parameter measurements 
since these assume the model is a perfect representation 
of the real galaxy. We list the GALFIT median measure- 
ment uncertainties in Table O Uncertainties for the low 
redshift cluster were slightly lower than for the two high 
redshift clusters. For all three clusters, uncertainties for 
r 1 / 4 law parameters were slightly lower than for the Sersic 
fits. However, this does not imply that these fits were nec- 
essarily better and could indicate that these models were 
a poorer representation of the real data. In cases where a 



model is a poor match to the data, statistical uncertainties 
can be smaller because a small ch ange in the parameters 
will produce a large change in \ 2 (Peng 20031) • The Ser- 
sic function with its extra free parameter is better able to 
fit a galaxy profile that deviates slightly from the tradi- 
tional r 1 / 4 law. While the uncertainties determined from 
profile fitting may be statistically accurate, they do not 
provide realistic uncertainties for the various derived pa- 
rameter values. We describe in the Section l5.1l our use of 
simulated galaxies to derive more realistic measurement 
uncertainties. 

4. RESULTS 

Our high redshift cluster galaxies have half-light radii 
for our assumed cosmology ranging in size from ~ 0.8 to 
35 kpc with magnitudes between 20.3 < i' < 23.7. The full 
RX J0152. 7-1357 spectroscopic sample includes 41 galax- 
ies, 36 of which are contained within the ACS images. 
Structural parameters have been obtained for all 29 spec- 
troscopically confirmed members. The RX J1226. 9+3332 
sample of 112 galaxies with measured parameters includes 
the original spectroscopic sample along with 63 galaxies 
having archival spectroscopy. 54 galaxies in this full sam- 
ple are members. The 28 galaxies in the RX J0142. 0+2131 
sample have magnitudes ranging from 16.9 < i' < 22.0 
and sizes 0.7 < r e < 34 kpc. In Tables [3] - [5j also avail- 
able electronically, we provide the measured and derived 
photometric and structural parameters for these galaxies. 
Galaxies are sorted by RA and columns are as follows: 
Column (1) Galaxy ID. This number comes from our orig- 
inal SExtractor detection in GMOS images. We use this 
number to identify the galaxies in all our publications. 
Column (2) RA (J2000.0) 
Column (3) DEC (J2000.0) 

Column (4) Number of measurements made from different 
images for each galaxy. Table entries are average values 
from all measurements. 

Column (5) Total apparent magnitude derived from r 1 / 4 
law profile fits. This is measured in the i'— band for RX 
J0152. 7-1357 and RX J0142.0+2131, and /-band for RX 
J1226. 9+3332. Magnitudes are uncorrected for reddening. 
Column (6) log (r e ) from our r 1 / 4 law fits, with r e in arc- 
sec. We take the radius as the geometric mean of the 
semi-major and minor axes, r = (a&) ' 5 . 
Column (7) Mean surface brightness within the effective 
radius derived from total magnitude and r e from our r 1 / 4 
law fits, in mag arcsec" 2 , using (fj,) e — m to t + 2.51og2-7r + 
2.5 log r-2. 

Column (8) Total magnitude from Sersic function fits, 
bands as in Column (5). 

Column (9) log (r e ) from Sersic function fits, with r e in 
arcsec. 

Column (10) (fi) e derived from total magnitude and r e 
from our Sersic function fits, in mag arcsec -2 . 
Column (11) Best fitting Sersic function parameter, n. 
Column (12) Position angle, North through East. 
Column (13) Ellipticity, derived from the fitted axial ratio. 
Column (14) Cluster members are denoted by a T', non- 
members by '0'. Not every galaxy in our sample proved to 
be a cluster member, but all sample galaxies for which we 
were able to measure structural parameters are included in 
these tables. A blank entry indicates no redshift is avail- 
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able. These include galaxy ID 910, required fitting as a 
neighbor of galaxy ID 899, galaxy IDs 1009 and 1253 which 
had archive spectroscopy but with too low S/N to derive 
a redshift, and galaxy ID 1254 for which we were unable 
to extract a redshift due to confusion with a second object 
in the slit. 

5. DISCUSSION 

5.1. Internal consistency 

Simulated galaxies are used to investigate sources of er- 
ror in our pipeline and test the accuracy of the software 
used. GALFIT, for example, produces uncertainties for 
the output parameters which tend to be too low since these 
are strictly random uncertainties and do not take into ac- 
count the fact that the profile we fit each galaxy with may 
not b e an exact match to the true galaxy profile (Peng 
2001). We therefore use the simulated galaxies to deter- 
mine more realistic measurement uncertainties and to test 
how these uncertainties affect our FP measurements. The 
simulated galaxies are also used to test the effect of dif- 
ferent PSF models and sizes for convolution in the galaxy 
profile fitting. In this section, we describe the methods 
used to simulate the galaxies, describe tests to determine 
the best PSF to use for the 2-D surface brightness profile 
fitting, and determine expected uncertainties in our galaxy 
parameter measurements. 

5.1.1. Galaxy simulations 

To generate realistic galaxies, we used the structural 
parameters of 148 galaxi es in our low redshift Com a clus- 
ter comparison sample ()j0rgensen fe Framd 11994? ). We 
transformed the values of Ms to the filter and red- 
shift of each clust e r usin g stellar population models of 
iBruzual fc Chariot! ((20 03?). For simulating galaxies at z 
= 0.83, we transform the Mb magnitudes of Coma cluster 
galaxies to the observed i' band using 

M B = B rest - DM(z), (6) 
B rest = i' + 0.8026-0.4268(i'-z')-0.0941(?;'-z') 2 , (7) 

(jJOrgensen et al.ll2005l l2007t ) and a mean empirical color, 
(i' — z') — 0.8, for our early type galaxies at this redshift. 
For z = 0.28, we use 

Brest =i' + 0.4753 +1.6421(r'-i')-0.0253(r'-i') 2 (8) 

(|Barr et alj|2005l ) with a mean empirical color (r' — i') = 
0.5. We do not include a separate simulation for z = 0.89 
as those made for z = 0.83 are considered representative. 
We scale i e ff to the appropriate redshift within our as- 
sumed cosmology. To create a large sample of simulated 
galaxies we added a small amount of random scatter to 
the magnitudes, ? e ff, and Coma galaxy axial ratios while 
randomly generating positions and position angles. The 
Sersic parameter, n, was allowed to vary randomly be- 
tween 0.5 - 5.5, except for a few sets of simulations where 
values ranged between 3.5 - 4.5 in order to increase the 
total number of early type galaxy simulations. The diski- 
ness/boxiness of each galaxy was also allowed to vary ran- 
domly. Galaxies were added as perfect Sersic profiles, al- 
beit with Poisson noise added to the galaxy images. 

In order to simulate galaxies with multiple components, 
we also generated a set of bulge + disk galaxies. In these 



cases, bulges were created by modeling from the Coma 
galaxy structural properties. Disks were then added to 
the same position assuming disk-to-total light ratios rang- 
ing randomly from 0.25 - 0.65 with disk sizes ranging from 
(1 - 1.4) r e (buige)- The inclination w as varied from to 
90 deg rees, uniformly in cosi. As in I Jorgensen fc Franxl 
(1994), we assign each component an intrinsic axial ratio: 
b/a(bulgc) = 0.7 and b/a(disk) = 0.15. These values arc 
chosen such that when transformed into observed values 
according to the inclination by 

b/a bserved = \A ~ (1 ~ (b / a intr i nsic ) 2 ) (cos i) 2 (9) 

(jSandage et al.|[l970[ ). they will display the same range of 
axial ratios found in real galaxies. 

Simulated galaxies were added to real images using two 
different means. Galaxies were generated with ARTDATA 
by defining a Sersic profile and using the Coma based pa- 
rameters. They were also created using GALFIT itself by 
turning off the fitting for all parameters thereby forcing 
GALFIT to output models of the input parameters. No 
differences were observed in the galaxies generated or in 
the results obtained from the two methods, but both were 
used to ensure that no bias was produced from the simu- 
lation method. 

Before being added to images, simulated galaxies were 
first convolved with a PSF and noise was added to the 
galaxies with the IRAF task MKNOISE. PSFs were cre- 
ated using 4-6 unsaturated real stars having high S/N us- 
ing the routines in the IRAF package DAOPHOT. These 
were created for each ACS image that simulated galaxies 
were added to. 

In order to measure realistic GALFIT parameter uncer- 
tainties we generated a total of over 2600 z = 0.83 and 
2500 z = 0.28 simulated galaxies. So as to not increase 
the surface density of objects in the images and thereby 
affect recovery of the surface brightness profile parameters, 
no more than 50 were added at a time to the real cluster 
images. It is important to note that all galaxies are gener- 
ated as perfect Sersic profiles and recovered as Sersic and 
r 1 / 4 law profiles. We do not simulate structure, such as 
arms, bars, double nuclei, shells, or twisted isophotes, all 
of which would effect the fitting. For the measurement 
of the cluster FPs, we have included only non-emission 
line, non-star forming galaxies which we take to be the 
early types. However, at higher redshift, one might expect 
to find more late-type features in these non-star forming 
galaxies. As a simple test of how added structure will af- 
fect single profile model fits, we use the set of simulated 
galaxies having both bulge + disk components. 

In Figure [4) we display a face-on view of the FP with 
our Coma cluster and high redshift galaxies plotted. We 
overlay our high redshift simulated galaxy sample. Ex- 
pected velocity dispersions given the input r e and (/z) e 
for each simulated galaxy were calculated using the FP 
equation for Coma, a = (log r e + 0.82 log( J) e + 0.443)/1.3 
(| Jorgensen et al.l l2006). It can be seen that our simulated 
galaxies span the same region of the FP as the real galax- 
ies. 

5.1.2. Choosing the best PSF for profile fitting 

In order to determine the most appropriate PSF to use 
for dcconvolution during the surface brightness profile fit- 
ting, we use a subset of 700 simulated galaxies to test the 
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recovery of parameters with GALFIT using several differ- 
ent PSF models and sizes. Images were run through our 
reduction pipeline starting with SExtractor object detec- 
tion in order to provide initial guesses for the GALFIT 
profile fitting. We tested the GALFIT recovery using four 
different PSFs: the PSF based on real stars, a raw 3" 
Tiny Tim PSF generated for the appropriate location on 
the ACS chip, and 3" and 9" drizzled PSFs (see Section 
13.31 for details) . Galaxies were fitted with both Sersic and 
r 1 / 4 law profiles. The r 1 / 4 law fits show greater discrep- 
ancy from input values than Sersic function fits because 
galaxies were added with a Sersic parameter n which var- 
ied randomly between 1 and 4.5 for these tests. 

In Table [5] we list the average differences and standard 
deviation between input and recovered magnitude, effec- 
tive radius (in arcsec), and the combination of parameters 
that enter into the FP, logr e — /31og(/) e , hereafter referred 
to as the Fundamental Plane parameter (FPP), where (I) e 
is the surf ace brightness within r e in units of Lq pc -2 and 

0.8 (jJorgensen et all [20061. The average difference 

between input and recovered values is nearly identical re- 
gardless of the PSF used. However, we do find that using 
the raw Tiny Tim PSF consistently provides the worst 
results. Using 9" drizzled PSFs showed very slight im- 
provement over 3" drizzled PSFs and very similar results 
to the real PSFs. A test using a subsampled PSF showed 
no improvement in the output results. Due to the dif- 
ficulties producing a real PSF from the few stars in our 
images, we therefore chose to use the model 9" drizzled 
PSF. In the event that the Tiny Tim modeled 9" drizzled 
PSF is not the best representation of the true PSF, we 
note that the differences in the FPP due to the PSF used 
are insignificantly small. 

5.1.3. Structural parameter measurement uncertainties 

The full set of simulated galaxies were used to investi- 
gate realistic uncertainties in structural parameter mea- 
surements. Average offsets (recovered - input values) and 
the associated rms scatter in these differences are provided 
in Tables [7] and [5] As with the PSF tests, simulated galax- 
ies were recovered using the same methods as real galaxies. 
The 9" drizzled PSF described in the previous section is 
used for the surface brightness profile fitting of both real 
and simulated galaxies. 

We compare the r 1 / 4 law recovered parameters from 
GALFIT with the input values for all simulated galaxies 
in Figure [5j We find small average systematic offsets and 
associated rms scatter from the true magnitude and log r e 
(arcsec) of —0.2 ± 0.4 and 0.1 ± 0.2 for the high z sample. 
For the lower redshift set, we find offsets of —0.4 ± 0.4 
and 0.3 ± 0.3 respectively. Average measurement offsets 
and rms scatter for the FPP are only —0.01 ± 0.06 and 
—0.04 ± 0.10 for the two samples. 

This is for the full set of simulated galaxies including 
those that were created with exponential profiles (n ~ 1) 
but which were recovered with n = 4, r 1 / 4 law profiles. 
To obtain realistic uncertainties in our FP parameters, we 
must compare only the range of parameters from galaxies 
that went into constructing our FPs. We plot a histogram 
of the Sersic n values measured by GALFIT for our real 
galaxy samples in Figure [5] The shaded histogram dis- 
plays only those galaxies used in our FP. Galaxies which 



were not included in our FP either had emission line spec - 
tra, log Mass < 10.3, or n < 1.5 (jJprgensen et al.ll2006l ). 
Wc therefore repeat the comparison between input and re- 
covered parameters in Figure [7J for galaxies created with 
n > 2.0 profiles. It can be seen from Table [7J that while 
magnitude and size errors have decreased, the systematic 
error and random uncertainty in the FPP of —0.01 ± 0.05 
and —0.02 ±0.10 for the two samples exhibit little change. 
The larger rms scatter for the lower redshift sample is due 
primarily to a few outliers with failed fits. In our real 
galaxy sample all failed fits are flagged and refitted or left 
out of the FP measurement. We therefore consider this 
measurement uncertainty to be an upper limit. 

Measurement errors are similar when fitting with a Ser- 
sic function, allowing n to vary. Systematic offsets are 
slightly lower for the magnitude and r e measurements, but 
there is no significant improvement in the FPP over fits 
made with r 1 ' 4 law profiles. In Figure [51 we display the 
error in recovered n as a function of input n. It is appar- 
ent that random error increases with increasing n. This is 
likely because fitting of n=l type galaxies tends to be more 
robust than n=4 galaxies where there is greater weight in 
the wings of the galaxy and fits are more easily affected 
by the sky values and near neighbors. 

Because we expect the intrinsic value of n to affect the 
results for r 1 / 4 law fits, we plot the error in r 1 / 4 law recov- 
ered parameters as a function of input Sersic n in Figure[5J 
Fluxes and sizes are overestimated for intrinsically low n 
galaxies because the r 1 ' 4 law profile imposes broader wings 
than the galaxy has. Magnitude and half-light radius are 
conversely underestimated at high n (> 4). Because flux 
and size are highly correlated, the errors in the FPP are 
reduced, but can be as large as 0.05 for n = 1 exponen- 
tial disk galaxies fitted incorrectly with a r 1 / 4 law profile. 
A linear fit to the offset from the true input values finds 
AFPP = 0.017 n - 0.067. 

We investigate other factors such as near neighbors 
which might affect the fits. Simulated galaxies with close 
neighbors within 1 r e or within 1 arcsec display a greater 
deviation in the recovered vs. input values for the FPP 
than those galaxies with a nearest neighbor at a distance 
greater than 2 arcsec. In Figure [TU1 we display the errors 
as a function of nearest neighbor distance. Points enclosed 
by triangles have neighbors within 1 r e . From the plot it 
is apparent that these points have a larger dispersion than 
those with more distant neighbors. The dispersion is 3 
times as large for objects having neighbors within 2 arcsec 
as compared to those with nearest neighbors > 6 arcsec. 
We also test whether errors change as a function of r e , 
magnitude, or surface brightness (Table [5]). We do find 
that random errors in the FPP may increase slightly at 
fainter magnitudes and smaller sizes. We find no trend 
as a function of the input c parameter, a measure of the 
galaxy core's intrinsic diskiness/boxiness. 

Finally, we investigate how single profile fits are affected 
by multiple component, bulge + disk, galaxies. Forcing a 
single r 1 / 4 law fit to a bulge embedded in a disk recovers 
a magnitude comparable to the sum of the two compo- 
nents and a half-light radius that is larger than for the 
bulge alone. When combined to form the FPP parameter, 
the average systematic offset from the FPP for the bulge 
component alone is 0.07 with a = 0.09 (Table [JJ. This 
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offset is larger than for the case of single profile galax- 
ies although similar to the random error in single profile 
cases. Obviously the extended exponential disk will affect 
the fit, but the resultant FPP is similar to that of the 
bulge FPP alone within the expected uncertainties. Thus, 
we expect a r 1 / 4 law fit to a distant galaxy with large bulge 
and low surface brightness disk will recover primarily the 
bulge component with only a slight offset in the FPP. Like- 
wise, we would expect the spectroscopic measurement of 
the velocity dispersion to come largely from light within 
the bulge. For cases where faint disks are not obvious, we 
expect FP measurements to primarily represent the bulge 
components. We conclude from our simulated galaxy anal- 
ysis that while r 1 / 4 law fits may lead to systematic errors 
in individual parameters for real galaxies, we do not ex- 
pect this to be a problem for our FP analysis. We explore 
this further in Section I5T21 

5.1.4. Comparing real galaxy measurements 

For both RX J0152. 7-1367 and RX J1226.9+3332, a 
number of galaxies were observed multiple times in over- 
lap regions from different telescope pointings. As one final 
test of our internal consistency we compare the derived 
structural and photometric properties from these multiple 
measurements. In Figure 1111 we compare the difference 
in derived FPP for both Sersic and r 1 ' 4 law profile fits in 
pairs of images. Different symbols denote different image 
pairs and the corresponding dashed lines show the aver- 
age difference for each pair. Average differences for pairs 
range from (absolute value) 0.0007 to 0.015 with standard 
deviations of 0.002 to 0.015. These differences are gen- 
erally smaller than GALFIT measurement uncertainties, 
and the scatter in these real measurement differences is 
smaller than the random uncertainties determined from 
simulated galaxies. We do not find any worrisome system- 
atic trends between images; all differences between images 
arc within expected measurement uncertainties. 

5.2. External consistency 

iBlakeslee et al.l (|2006l ) have used the same ACS dataset 
as used in this paper to study the structural parameters 
and colors of 149 galaxies in RX J0152. 7-1357. They also 
perform surface brightness profile fitting using GALFIT 
with a procedure very similar to ours. Images are pro- 
cessed in a similar manner, and like us, they fit all neigh- 
bors with ipnsw < 25 AB while masking objects fainter 
than this. Differences in methodology are few but three 
significant differences exist. While we use a 9" Tiny Tim 
model PSF drizzled in the same manner as our data, they 
choose to use empirical PSFs from archival HST images. 
They also impose an upper limit for the fitted Sersic pa- 
rameter, n, constraining the value to n < 4, whereas we 
applied no constraints. Because of the strong coupling 
between Sersic function parameters, this will cause differ- 
ences in the measured values of both r e and total magni- 
tude for galaxies we find with n > 4. Finally, while we hold 
the sky value constant during the fitting at the SExtractor 
measured local value, they simultaneously fit for the sky 
value with GALFIT, except in cases where this led to poor 
fits. 

We find a total of 26 objects in common, 3 of which 
have large discrepancies in magnitude. These were galax- 
ies that were not well fit with single Sersic functions and 



which had late type features evident in the residual im- 
ages. Excluding these, we compare the measured Sersic 
surface brightness profile parameters of the remaining 23 
galaxies to the values listed in the Blakeslee et al. cat- 
alog. We plot the differences in our measurements in a 
Alogr e vs. A(//) e plane (Figure [T2"|) . From a linear fit 
to these data points we measure a magnitude zeropoint 
shift of —0.025 between the two works from the offset in 
the intercept from zero. This is in the sense that we find 
objects to be on average 0.025 mag brighter. We find that 
large differences in measured surface brightness or size are 
correlated with large discrepancies in the recovered Sersic 
n parameter (see also Figure [5]). Galaxies with measured 
n differing by less than 0.1 for the two groups are circled. 
Correcting for the zeropoint shift we compare the differ- 
ence in magnitudes, r e , n, and the FPP. Results are listed 
in Tableland displayed in Figure 

The increasing difference in measured n at n > 4 is eas- 
ily explained by the constraints Blakeslee et al. impose 
on the Sersic n. Blakeslee et al. chose not to use Tiny 
Tim models finding that analytical representations for the 
PSF overestimated fitted r e values. However, with our 9" 
model PSF we find r e values 2% smaller on average than 
their measurements. Average differences in total measured 
magnitude and effective radius are insignificantly small 
compared to the rms scatter. However, there appears to 
be a trend as a function of size and magnitude. We find 
brighter galaxies to be brighter than found by Blakeslee et 
al. and fainter galaxies to be fainter. Likewise, we measure 
larger sizes for larger galaxies and smaller sizes for smaller 
galaxies in comparison to Blakeslee et al. This can be at- 
tributed in part to different values recovered for the Sersic 
n parameter. In five cases where we find similar values 
within 0.1 for measured n, magnitude and size discrepan- 
cies are negligible regardless of size or magnitude. For 10 
objects measured by Blakeslee et al. to have n = 4.0 due 
to their constraint on the Sersic index, and for which we 
find best Sersic function fits with n > 4.0, we compare the 
Blakeslee et al. results to our r 1 / 4 law (n — 4.0) fits. We 
find the average offsets in magnitude and size are similar 
but the random scatter is nearly half that of the differences 
in Sersic measurements for these same galaxies. For large 
galaxies, in cases where we find n > 4.0, the covariance 
of Sersic function parameters generates fits with larger r e 
than found by Blakeslee et al. For small galaxies which are 
more strongly affected by the choice of PSF, we speculate 
that our broad 9" PSF may be the cause of the smaller 
size measurements as compared to Blakeslee et al. Due 
to the coupling of Sersic function parameters, a smaller 
recovered size results in a fainter recovered magnitude. In 
Figure [T3] we also show the parameter measurement errors 
in the A(/i) e — Alogr e plane for our z = 0.83 simulations 
fitted with an r 1 / 4 law profile. We find the same coupling 
of errors and an offset in magnitude of 0.018, similar to 
that found for the Blakeslee et al. comparison. Since we 
used PSFs generated from real stars to create the simu- 
lated galaxies and the 9" model PSF to fit the galaxies, 
this further suggests that the broad PSF may be cause of 
the slight magnitude differences. 

Because of the correlation between Sersic profile param- 
eters, the FPP is insensitive to these small variations in 
size measurement. When taking into account the small 
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magnitude zero point differences, we find the differences 
between the Blakeslee et al. and our FPP measurement to 
be completely negligible and smaller than the systematic 
error of 0.006 determined from galaxy simulations. The 
rms scatter in the differences between our measurements 
is a factor of 5 smaller than the uncertainty expected from 
the simulations. The systematic difference in the FPP 
measurement of 0.009 is larger when not including the zero 
point shift, but is still much smaller than our expected ran- 
dom uncertainty of 0.05 determined from simulations. 

5.3. Implications for FP measurements 

Although we have fit each galaxy with both Sersic and 
r 1 / 4 law functions, we use measurements of r e and surface 
brightness from the r 1 / 4 law fits to derive the FP since 
the galaxies in our low redshift comparison sample were 
originally fitted only with r 1 / 4 law profiles. Not all of our 
sample consists of pure r 1 / 4 law galaxies, and the Sersic 
function fits find a wide range of values for the Sersic pa- 
rameter. Howeve r, only galaxies with n > 2 are used in 
our F P analysis (|j0rgensen et alj 120061 l2007t iBarr et alj 
2006j). In this section we explore how galaxy parameter 
measurement uncertainties translate to FP measurement 
uncertainties and whether any systematic biases exist in 
our methodology that would affect the derivation of the 
FP. 

A plot of magnitude and FPP differences between GAL- 
FIT Sersic and r 1 / 4 law fits of real galaxies vs Sersic in- 
dex, n, shows very good agreement for n ~ 4 galaxies as 
expected (Figure IT4"|) . Many galaxies were best fit with n 
< 4 and, for these, the r 1 / 4 law measured magnitude was 
as much as 1 mag brighter than that derived from the Ser- 
sic fits. Although output parameters from Sersic and r 1 / 4 
law fits often differed, a plot of the difference in magnitude 
and in surface brightness vs. difference in log r e for the 
two fits shows the parameters to be strongly correlated and 
explains why differences in the FPP, although also corre- 
lated with Sersic n, are smaller. Galaxies with best Sersic 
fit n > 2.0 are displayed as larger points, but note that 
even galaxies with low Sersic index follow the same linear 
relation. Indeed, the FPP is reasonably robust against the 
form of the pr ofile used to fi t the g a laxies. This ef f ect ha s 
bee n noted by iKelson et all (|2000li . iTruiillo et all (|200lD . 
and iLucevi (|1997l ). While the uncertainties in measured 
structural and photometric parameters may seem large in 
cases where n < 4.0 according to best Sersic function fits, 
the FPP has a much lower ~ 0.01 systematic error with 
0.1 random uncertainty. We find that even when we fit a 
pure exponential disk (n = 1) with an r 1 / 4 law profile, the 
errors in the magnitude and r e are large but the average 
error in the FPP is only 0.05. 

In our study of the FP we investigate whether the slope 
changes as a function of redshift. A change in slope would 
imply a change in M/L ratio as a function of galaxy mass. 
It is therefore critical to ensure our measurements arc not 
biased as a function of magnitude or size. From our galaxy 
simulations we find some evidence for a slight increase of 
up to a factor of 3 in the scatter of the FPP at smaller 
sizes (< 1 arcsec) and at fainter magnitudes (i > 22.5). 
However, this effect is small. More significantly, wc find 
a strong correlation between the uncertainty in the FPP 
and the difference in measured vs. intrinsic Sersic index 



n. Exponential disk galaxies with n < 2 fitted with r 1 / 4 
law profiles have systematic offsets as large as the random 
uncertainty. If there is a trend in real galaxy magnitudes 
with index n, this could affect the FP slope measurement. 
In Figure [TBI we display the measured Sersic index n ver- 
sus SExtractor total magnitude for the galaxies in our RX 
J0152. 7-1367 and RX J1226.9+3332 FP samples. For the 
most part, the points exhibit random scatter. There is a 
hint of a trend in that the faintest galaxies have a slightly 
lower average Sersic index n. The correlation coefficients 
for these two samples are small, only -0.31 and -0.14 for 
RX J0152. 7-1367 and RX J1226.9+3332 respectively. Al- 
though this slight correlation coupled with the trend of in- 
creasing FPP measurement offset with decreasing n may 
affect the FP slope, this effect is much smaller than the 
random uncertainties in the FPP. 

It is also important to ensure our measurement uncer- 
tainties are insensitive to redshift. A redshift bias could 
be caused by, for example, GALFIT parameter measure- 
ment uncertainties being systematically different for ob- 
jects with larger sizes or brighter magnitudes at low red- 
shift as compared with smaller, fainter objects at high red- 
shift. We do not find any evidence for this in our recovery 
of simulated galaxy parameters in simulations of the high 
(z = 0.83) and low (z = 0.28) redshift clusters. The average 
magnitude and effective radius measurement uncertainties 
are slightly larger for the low redshift sample but scatter in 
the recovered values is about the same for both clusters. 
The FPP has an average offset of recovered values from 
input which is insignificant for both populations, and a 
random uncertainty which is larger for the lower redshift 
cluster due to a small number of outliers with particu- 
larly poor fits. In our real galaxy sample all failed fits are 
flagged and refitted or left out of the FP measurement. 

Through this simulated galaxy analysis we have found a 
number of other factors which affect the galaxy structural 
parameter measurements. Near neighbors affect profile fits 
and based on simulations, we find that the highest surface 
brightness, most compact galaxies are affected the least. 
However, the only effect this has on lower surface bright- 
ness galaxies is to increase the random uncertainty by up 
to factor of 4 for only those few galaxies which have a 
neighbor within 1 arcsec or 1 r e . No systematic variation 
as a function of surface brightness is expected. 

Although we cannot simulate galaxies with a full range 
of substructures, nor do we know the true distribution of 
such features in our real data, we estimate uncertainties 
for galaxies with more complex structures with bulge + 
disk simulations. We find a systematic offset of 0.07, much 
larger than the 0.01 for simulated galaxies generated with 
a single profile, with a slightly larger rms scatter of 0.10. 
With a large sample of galaxies exhibiting complex mor- 
phologies, we would expect measurements of the FP to be 
affected by the larger uncertainties. However, we check 
the goodness of all fits from residual images. Any fea- 
tures not modeled by the simple Sersic/r 1 / 4 law profile fits 
will show up in these images as imperfect subtraction from 
the galaxy model. For most of our galaxies, subtraction 
is good without any obvious sign of extra arms, twisted 
isophotes, low surface brightness disks, or any other struc- 
tures. Galaxies with such features tend to be best fit with 
Sersic n < 2.0, and these objects are not included in our 
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FP measurements. 

In conclusion, we take final uncertainties in the FPP to 
be ~ 0.1 based on our simulations, dominated by random 
uncertainties. We find no systematic offsets between re- 
covered and input values for the FPP as a function of 
magnitude or half-light radius from our galaxy simula- 
tions. While we do find uncertainties are affected by near 
neighbors or Sersic index n, we do not expect the shape 
or tilt of the FP to be affected by these uncertainties, un- 
less for example all lower mass galaxies have exponential 
profiles, which we see no evidence for. We find that our 
measurement uncertainties should not vary with redshift 
unless, for example, higher redshift galaxies contain more 
late-type morphological features. Although this is a possi- 
bility, we would expect to spot this in our residual fitting 
images. We therefore do not find any evidence for mea- 
surement uncertainties or systematic biases introduced by 
our data reduction methods to affect our derived FP as a 
function of redshift, galaxy size, or magnitude. We expect 
the measured tilts to be intrinsic values for each cluster. 

6. SUMMARY 

Using new and archival HST/ACS high resolution im- 
ages observed in the F775W or F81AW bands, we have 
measured structural and photometric properties for galax- 
ies in the clusters RX J0142.0+2131, RX J0152. 7-1367, 
and RX J1226. 9+3332. These data are being used to in- 
vestigate the Fundamental Plane as a function of redshift, 
mass, and environment as part of our Gemini/HST Galaxy 
Cluster Project aimed at understanding cluster galaxy evo- 
lution. In this work we have described our data processing 
and analysis methods, presented catalogs of parameters 
measured from 2-D surface brightness profile fitting, and 
discussed expected uncertainties in our measurements and 
in the quantities which go into the measurement of the FP. 

We determine average uncertainties in measured param- 
eters for n > 2 galaxies from galaxy simulation tests. We 
define systematic errors and random uncertainties as the 
mean and rms scatter in the difference between recovered 
and input simulation values. For measured total magni- 
tudes, these are —0.1 ± 0.3 for r 1 / 4 law profiles. The aver- 
age error in size, log r e , is 0.1 ±0.2. Random uncertainties 
are larger than any systematic offsets from intrinsic values. 
Sersic profile fits do slightly better than r 1 ' 4 law profiles in 



recovering the input values of the simulated galaxies with 
magnitude errors of ~ 0.03 ± 0.3 and errors in log r e of 
~ —0.02 ± 0.2. However, the rms scatter is similar for fits 
with either profile. 

Due to the combination of parameters that enter the FP, 
log r e + /31og(7) e , we find structural and photometric un- 
certainties to have little effect on the FP measurement. We 
take final uncertainties in this FP parameter to be ~ 0.1, 
dominated by random error. We find no systematic offsets 
between recovered and input values for the FPP in our 
galaxy simulations as a function of magnitude or effective 
radius. While we do find uncertainties are affected by near 
neighbors or Sersic index n, we do not expect the shape 
or tilt of the FP to be affected by these uncertainties, un- 
less for example all lower mass galaxies have exponential 
profiles, something not seen in our samples. We find that 
our measurement uncertainties should be invariant with 
redshift unless, for example, higher redshift galaxies con- 
tain more late-type morphological features. While this is 
a possibility, we would expect to spot this in our residual 
fitting images, so do not expect this to significantly affect 
our results. We conclude that the methods used to derive 
the effective parameters r e and (fi) e do not introduce any 
significant bias in our FP measure ments, or in the r esults 
prese n ted for thes e three clusters in lJorgensen et all ([2006, 
l2007h : lBarr et all (|2006h . 
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Table 1 







HST/ACS Imaging 


Data 




Cluster 


# fields filter 


total texp(s) 


Program ID 


RXJ0142. 0+2131 


2 


F775W (i') 


4420 


9770 


RXJ0152. 7+1357 


4 


F625W (r') 


4750 


9290 




4 


F775W (i') 


4800 






4 


F850LP (z') 


4750 




RXJ1226.9+3332 


4 


F606W (i?) 


4000 


9033 




4 


F814W (J) 


4000 
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Table 2 



GALFIT Median Measurement Uncertainties 





z= 0.83-0.89 




z= 0.28 






Sersic fit 


r l/4 


Sersic 


r l/4 


Magnitude 


0.010 


0.006 


0.003 


0.001 


r e 


0.009 


0.009 


0.009 


0.008 


(M>e 


0.040 


0.025 


0.013 


0.008 



Note. — Uncertainties are similar for the two high red- 
shift clusters and arc therefore combined in this table. 



Table 3 



RXJ0152. 7-1357: Photometric and Structural Parameters 



ID 


RA 


Doc 


N meas 


Hltot,dev 


log(r e )deu 






log(r c )ser 




n 3e r 


PA„„ e 




Mcmber a 




(2000) 






(F775W) 


(arcscc) 


mag arcsec - 2 


(F775W) 


(arcscc) 


mag arcsec - 2 










193+ 


1:52:50.80 


-13:55:28.9 


2 


20.25 


0.115 


22.82 


20.64 


-0.161 


21.83 


2.13 


34.3 


0.07 





264 


1:52:44.66 


-13:55:37.3 


1 


21.17 


-0.181 


22.26 


21.17 


-0.167 


22.89 


4.33 


-29.2 


0.47 





338 


1:52:43.33 


-13:55:44.4 


1 


21.93 


-0.447 


21.69 


22.19 


-0.632 


21.02 


2.18 


-17.4 


0.33 


1 


346 


1:52:37.42 


-13:55:50.1 


2 


21.11 


-0.408 


21.06 


21.21 


-0.480 


20.80 


3.29 


133.8 


0.25 


1 


422 


1:52:34.59 


-13:55:58.8 


1 


21.80 


-0.494 


21.32 


21.79 


-0.487 


21.34 


4.06 


73.7 


0.09 


1 


460 


1:52:36.11 


-13:56:08.5 


1 


20.87 


-0.365 


21.04 


20.91 


-0.387 


20.96 


4.16 


101.9 


0.25 





523 


1:52:42.38 


-13:56:18.7 


2 


21.10 


-0.360 


21.30 


21.20 


-0.415 


21.12 


4.48 


-10.2 


0.14 


1 


566 


1:52:38.03 


-13:56:28.1 


2 


20.90 


-0.097 


22.41 


20.75 


0.010 


22.79 


4.54 


1.5 


0.53 


1 


627 


1:52:38.48 


-13:56:33.6 


2 


21.81 


-0.350 


22.05 


21.88 


-0.388 


21.93 


4.18 


-2.3 


0.25 


1 


643 


1:52:45.60 


-13:56:40.0 


2 


21.41 


-0.168 


22.58 


21.48 


-0.218 


22.38 


3.58 


132.5 


0.46 


1 


737 


1:52:45.77 


-13:56:46.1 


2 


21.85 


-0.329 


22.20 


21.90 


-0.365 


22.06 


3.71 


54.2 


0.08 


1 


766 


1:52:45.83 


-13:56:59.2 


2 


20.33 


0.063 


22.63 


20.33 


0.057 


22.61 


3.96 


32.6 


0.22 


1 


I IV) 


1 . ro. Oft A ft 


1 o . re. eg c 
-lo.OD.Oi.O 




21 44 


- U . o 1 D 


1 Q ir. 


21 45 


n fi9 ^ 

- U . O A O 


1Q OO 

ly.oo 


O.OD 


r.A ft 


U.Dl 




813 


1 -59-44 Q7 


1 3-57-04 9 

J-O . O l . IJ-±. _ 


4 


20.68 


-0 209 


21 63 


20 49 


-0 051 


22 22 


5 20 


5 


64 


\ 


896+ 


1:52:3b. 99 


i o.ct in i 

-13:57:10.1 


1 


22.29 


-0.179 


23.39 


21.87 


-0.234 


22.68 


0.31 


-13.3 


0.32 


1 


908 


1:52:43.74 


-13:57:19.4 


4 


20.80 


0.074 


23.16 


20.78 


0.067 


23.11 


3.80 


50.9 


0.21 


1 


1027 


1:52:43.32 


-13:57:2b. 7 


4 


21.83 


-0.413 


21.76 


21.94 


-0.489 


21.49 


3.40 


13.3 


0.24 


1 


1085 


1:52:42.94 


-13:57:35. U 


4 


21.23 


-0.273 


21.86 


21.15 


-0.208 


22.10 


4.52 


72.2 


0.27 


1 


1110 


1:52:39.93 


-13:57:42.6 


2 


21.09 


-0.143 


22.37 


21.23 


-0.254 


21.95 


3.15 


-26.2 


0.26 


1 


1159+ 


1:52:36.18 


-13:57:48.8 


1 


21.61 


-0.411 


21.55 


21.73 


-0.481 


21.32 


4.45 


-7.9 


0.47 


1 


1210 


1:52:42.83 


-13:57:55.3 


2 


21.81 


-0.390 


21.86 


21.98 


-0.508 


21.43 


2.87 


122.5 


0.29 


1 


1245 


1:52:43.57 


-13:58:00.0 


1 


21.18 


0.064 


23.50 


21.93 


-0.416 


21.84 


0.88 


21.6 


0.27 





1299+ 


1:52:47.34 


-13:59:26.1 


1 


21.28 


-0.071 


22.92 


21.28 


-0.134 


22.60 


2.49 


-22.3 


0.48 


1 


1385+ 


1:52:39.36 


-13:59:04.5 


1 


22.10 


0.021 


24.20 


21.69 


-0.209 


22.64 


1.11 


79.4 


0.34 


1 


1458 


1:52:39.64 


-13:58:56.6 


2 


21.94 


-0.541 


21.22 


21.97 


-0.547 


21.23 


4.47 


21.2 


0.29 


1 


1494 


1:52:39.08 


-13:58:48.8 


2 


18.44 


-0.252 


19.17 


18.41 


-0.223 


19.28 


4.44 


108.0 


0.51 





1507 


1:52:34.48 


-13:58:42.2 


1 


22.14 


-0.436 


21.94 


22.17 


-0.460 


21.85 


3.76 


61.6 


0.56 


1 


1567 


1:52:39.62 


-13:58:26.7 


2 


20.17 


0.472 


24.51 


20.43 


0.292 


23.88 


3.17 


50.4 


0.20 


1 


1590 


1:52:38.87 


-13:58:32.0 


2 


22.09 


-0.531 


21.42 


22.22 


-0.629 


21.06 


2.98 


85.7 


0.34 


1 


1614 


1:52:51.96 


-13:58:17.1 


1 


20.88 


0.036 


23.05 


20.80 


0.098 


23.28 


4.43 


62.1 


0.15 


1 


1682 


1:52:51.96 


-13:58:15.6 


1 


21.16 


-0.163 


22.34 


21.54 


-0.446 


21.30 


1.64 


127.7 


0.49 


1 


1811 


1:52:38.63 


-13:59:20.8 


1 


22.43 


-0.495 


21.94 


22.63 


-0.644 


21.40 


2.62 


3.3 


0.37 


1 


1920+ 


1:52:39.70 


-13:59:14.3 


1 


21.94 


-0.013 


23.87 


21.84 


-0.160 


23.03 


0.85 


-33.4 


0.92 


1 


1935 


1:52:41.88 


-13:59:53.6 


1 


21.36 


-0.295 


21.87 


21.29 


-0.242 


22.07 


4.57 


10.4 


0.18 


1 


1970+ 


1:52:48.03 


-13:59:58.6 


1 


19.78 


0.389 


23.72 


20.58 


-0.166 


21.74 


0.62 


57.4 


0.78 





2042+ 


1:52:42.38 


-13:59:46.6 


1 


18.12 


0.958 


24.90 


19.43 


0.212 


22.48 


0.75 


-40.3 


0.17 






a l: galaxy is considered a member of RX J0152. 7-1357; 0: galaxy is not a member 
"^Galaxy exhibits late-type structure and is not well fit with an r 1 / 4 law profile. 
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a l: galaxy is considered a member of RX J1226. 9+3332; 0: galaxy is not a member, left blank indicates rcdshift not available 
* Galaxy exhibits late-type structure and is not well fit with an r 1//4 law profile. 



Table 5 



RXJ0142. 0+2131: Photometric and Structural Parameters 
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\M/ e t dev 

mag arcscc 


-2 


mtot ser 


1 n < >' { r 1 

Iu & V ejser 

(arcscc) 


(M) e,ser 

mag arcsec" 


-2 




PA ave 


a ' 


ue 


]VI c nib c r 


1 


1:42:09.11 


21 


:33:23 


84 


1 


17.46 


0.360 


21 


.24 


16.93 


0.794 




22 


89 


6.98 


-57.7 





.07 


i 


22 


1:42:08.68 


21 


:33:22 


.62 


1 


18.98 


-0.313 


19 


.40 


18.97 


-0.305 




19 


.43 


4.09 


-85.9 





.28 


i 


88 


1:42:09.21 


21 


:33:13. 


.97 


1 


21.74 


-0.611 


20 


.67 


22.04 


-0.814 




19 


.96 


1.48 


53.8 





.50 


i 


128 


1:42:07.37 


21 


:33:02 


.13 


1 


18.76 


0.029 


20 


.90 


18.80 


0.002 




20 


.80 


3.77 


-0.2 





.21 


i 


205 


1:42:04.40 


21 


:32:39 


.82 


1 


18.47 


0.015 


20 


.54 


18.24 


0.209 




21 


.28 


5.72 


39.4 





.34 


i 


318 


1:41:57.29 


21 


:32:27 


.84 


1 


19.73 


-0.058 


21 


.43 


20.04 


-0.294 




20 


.56 


1.91 


-14.7 





.41 


i 


322 


1:42:03.69 


21 


:32:15 


.34 


1 


18.58 


0.070 


20 


.92 


18.85 


-0.142 




20 


.14 


1.91 


-19.2 





.65 


i 


379 


1:42:01.92 


21 


:32:10 


.19 


1 


18.71 


-0.276 


19 


.32 


18.68 


-0.252 




19 


.41 


4.30 


-37.4 





.52 


i 


412 


1:42:02.45 


21 


:31:57 


.59 


1 


18.25 


-0.061 


19 


.93 


18.19 


-0.019 




20 


.09 


4.40 


-47.0 





09 


i 


442 


1:42:02.64 


21 


:32:06 


.54 


1 


20.74 


-0.539 


20 


.04 


20.87 


-0.634 




19 


.69 


2.77 


-33.0 





.32 


i 


479 


1:42:03.46 


21 


:31:17. 


.36 


2 


16.31 


0.925 


22 


.92 


16.32 


0.919 




22 


.91 


4.01 


-58.6 





.38 


i 


537 


1:42:08.64 


21 


:31:45 


.50 


1 


19.97 


-0.461 


19 


.65 


19.86 


-0.370 




20 


.00 


5.07 


6.6 





.21 


i 


614 


1:42:01.26 


21 


:31:31 


.98 


1 


19.68 


-0.525 


19 


.05 


19.77 


-0.594 




18 


.80 


2.95 


52.3 





.33 


i 


637 


1:42:01.38 


21 


:31:22 


.13 


1 


19.41 


-0.246 


20 


.17 


19.43 


-0.266 




20 


.09 


3.78 


19.5 





.32 


i 


671 


1:42:03.21 


21 


:31:11 


.90 


2 


17.90 


-0.142 


19 


.18 


18.00 


-0.223 




18 


.88 


3.24 


25.7 





.60 


i 


760 


1:42:01.28 


21 


:31:04 


.59 


2 


19.87 


-0.566 


19 


.03 


19.93 


-0.608 




18 


.88 


3.22 


3.7 





.48 


i 


777 


1:41:59.76 


21 


:30:57 


91 


2 


19.11 


-0.364 


19 


.28 


19.16 


-0.401 




19 


.14 


3.54 


-6.9 





.44 


i 


844 


1:42:07.17 


21 


:30:49 


.71 


1 


19.69 


-0.437 


19 


.49 


19.65 


-0.409 




19 


.60 


4.37 


-76.7 





.47 


i 


911 


1:42:03.11 


21 


:30:31 


.59 


1 


20.67 


-0.579 


19 


.76 


20.65 


-0.561 




19 


.83 


4.25 


80.1 





.31 


i 


1012 


1:42:01.75 


21 


:30:17 


.39 


1 


20.23 


-0.500 


19 


.72 


20.19 


-0.473 




19 


.82 


4.35 


73.9 





.13 


i 


1029 


1:41:55.20 


21 


:30:12 


.08 


1 


19.85 


-0.278 


20 


.46 


19.82 


-0.249 




20 


.56 


4.25 


45.0 





.25 


i 


1043 


1:41:58.57 


21 


:30:01 


91 


1 


18.69 


-0.153 


19 


.91 


18.70 


-0.163 




19 


.87 


3.89 


5.7 





.50 


i 


1099 


1:42:05.63 


21 


:30:03 


.39 


1 


21.06 


-0.068 


22 


.71 


21.58 


-0.451 




21 


31 


1.06 


52.2 





.73 


i 


1179 


1:42:00.91 


21 


:29:41 


.60 


1 


18.93 


0.352 


22 


.68 


19.43 


-0.014 




21 


.35 


1.37 


-26.2 





.29 


i 


1205 


1:41:53.41 


21 


:29:26 


.73 


1 


19.34 


-0.351 


19 


.57 


19.33 


-0.342 




19 


.60 


4.10 


-19.3 





.64 


i 


1207 


1:42:04.04 


21 


:29:35 


.51 


1 


19.81 


-0.403 


19 


.78 


19.77 


-0.373 




19 


.90 


4.33 


60.2 





.46 


i 


1412 


1:42:07.28 


21 


:28:56 


.51 


1 


19.79 


-0.239 


20 


.59 


19.64 


-0.138 




20 


.95 


4.80 


25.6 





.14 


i 


1416 


1:42:06.40 


21 


:28:38 




1 


19.99 


0.155 




.76 


20.43 


-0.194 




21 


.45 


1.61 


-26.8 





.70 


i 



a l: galaxy is considered a member of RX J0142. 0+2131; 0: galaxy is not a member 
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Table 6 

PSF Tests (Recovered - Input values) 



Parameter 


Real a 


PSF 
Raw TinyTim 


drizzled 3" 


drizzled 9" 


(A(logr e + 0.81og(7) e )) rl/ 4 

(A(logr e + 0.81og(/) e )) ser 

(Ama#) rl/4 

(Amag) ser 

(Alogr e ) rl/4 

(Alogr e ) ser 


0.03 ±0.09 
0.01 ±0.05 

0.2 ±0.3 
-0.1 ±0.2 
-0.2 ±0.2 

0.0 ±0.1 


0.05 ±0.10 
0.01 ±0.04 

0.3 ±0.3 
-0.1 ±0.1 
-0.2 ±0.2 

0.1 ±0.1 


0.03 ±0.10 
0.01 ±0.06 

0.2 ±0.3 
-0.1 ±0.2 
-0.1 ±0.2 

0.0 ±0.1 


0.02 ±0.09 
-0.01 ±0.04 

0.2 ±0.3 
-0.0 ±0.2 
-0.1 ±0.2 

0.0 ±0.1 



Note. — Values listed are the average offset and standard deviation of the measurements. 
a Empirical PSFs generated from real stars 
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Table 7 

Profile Fitting Tests (Recovered - Input values) 



Parameter 


z=0.83 
All sims 


z=0.28 
All sims 


z=0.83 

User ^* 2 


z=0.28 

^ser ^* 2 


Amag r i/4 
Alogr e rl/4 

A(/i)e rl/4 

A FPP rl/4 


-0.19 ±0.41 
0.12 ±0.24 
0.43 ±0.85 

-0.01 ±0.06 


-0.38 ±0.43 
0.25 ±0.34 
0.87 ± 1.41 

-0.04 ±0.10 


-0.08 ±0.27 
0.05 ±0.17 
0.18 ±0.62 

-0.01 ±0.05 


-0.12 ±0.28 
0.09 ±0.21 
-0.32 ± 1.02 
-0.02 ±0.10 


Amag ser 
A log r e ser 

A(/i) e ser 

A FPP ser 


0.01 ±0.30 
-0.01 ±0.19 
-0.06 ±0.73 
0.00 ±0.05 
-0.2 ± 1.1 


-0.04 ±0.46 
0.02 ±0.23 
0.03 ± 1.10 
0.00 ±0.09 
0.2 ± 1.7 


0.02 ±0.32 
-0.01 ±0.21 
-0.07 ±0.81 
0.01 ±0.05 
-0.2 ± 1.3 


-0.04 ±0.52 
0.02 ±0.23 
0.03 ± 1.20 
0.01 ±0.10 
0.2 ± 1.9 


AmagBDrl/4 
Al0gr e BDrl/i 
A(M)e S_Drl/4 

A FPP BI>1 / 4 


-0.56±0.67 
0.19 ±0.24 
0.37 ±0.61 
0.07 ±0.09 









Note. — rl/4 refers to profile fits with r 1 / 4 law profiles, ser to Sersic profile 
fits, and BDrl/4 to objects generated with 2 components but fit with a single 
r 1 / 4 law profile. Errors listed are the rms scatter in the differences. 
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Table 8 

r 1/4 Law Profile Fitting Tests (Recovered - Input values) 





range 1 


(AFPP) 
range2 


range3 


Nearest neighbor 


d < 2" 


2 < d < 6" 


d > 6" 




-0.01 ±0.12 


-0.00 ±0.03 


0.00 ±0.04 


Size 


logr e > 


-1 < logr e < 


logr e < — 1 




-0.00 ±0.05 


-0.01 ±0.05 


-0.01 ±0.10 


Mag 


i' < 21 


21< i' < 22.5 


i' > 22.5 




0.00 ±0.04 


-0.00 ±0.03 


-0.01 ±0.07 


We 


<ji) e < 20 


20 < (fi) e < 22.5 


(li)e > 22.5 




-0.00 ±0.08 


-0.01 ±0.05 


-0.00 ±0.05 
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Table 9 

Profile Fitting Tests (external comparison) 



Parameter 


(This work - Blakeslee et al.) a 




(subset Blakeslee et al. 


n = 4) b 




zero point shift applied 


no shift 


This work: Sersic 


r 1 / 4 fits 


Amag ser 


0.01 ±0.08 


-0.02 ±0.08 


-0.01 ±0.09 


0.01 ±0.05 


A log r e ser 


-0.01 ±0.06 


-0.01 ±0.06 


0.01 ±0.07 


-0.02 ±0.03 


^(/^)e ser 


-0.02 ±0.21 


-0.04 ±0.21 


0.04 ±0.26 


-0.07 ±0.10 


An 


0.2 ±0.4 


0.2 ±0.4 






A FPP ser 


0.00 ±0.01 


0.01 ±0.01 


-0.003 ±0.014 


0.006 ± 0.006 



a Comparison between 23 galaxies in common. 

b Comparison between 10 galaxies for which we find best Sersic fits with n > 4.0 while Blakeslee et al. 
limit n < 4.0. We show the comparison with both our Sersic and r 1 / 4 law (n = 4.0) fits. 
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Fig. 1. — Several examples of output images produced by GALFIT. These include, from top to bottom, the panel from the original image, 
a model image with the primary galaxy fit with an r 1 ' 4 profile, and the residuals from the best fit. In every case, the primary galaxy being 
fitted is the central object. Galaxies are from RX J0152.7-1357: 737 (with Scrsic fit n = 3.71), 643 (n = 3.58), 1567 (n = 3.17), 1385 (n = 
1.11). Image sections are 5.2 arcsec across. 




Fig. 2.— Best r 1 / 4 law fits as in Figured] Galaxies are from RX J1226.9+3332. From left to right: 
= 4.74), and 563 (n = 3.92). Image sections are 11.8, 5.9, and 5.9 arcsec across respectively. 



1025 (with Sersic fit n = 6.96), 529 (n 
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Fig. 3.— Best r 1 / 4 law fits as in Figure [T] A further 3 examples from RX J1226.9+3332: 630 (with Sersic fit n 
and 1251 (n = 1.90). Image sections are 11.8 arcsec across. 



= 3.29), 739 (n = 3.14), 
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-0.5 0.5 1 1.5 

(2.211og r. - O.B210g<[>, + 1.841og(o)) / 2.66 



Fig. 4. — The face-on FP. Blue squares represent galaxies of RX J0152. 7-1357, green squares RX J1226. 9+3332, and red dots Coma sample 
galaxies. Black stars denote the location of our simulated galaxies in this FP. These simulated galaxies span the full range of the high z FP. 
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Fig. 5. — Left: z = 0.83 simulated galaxy set. Right: z = 0.28 set. From fits with r 1 / 4 law profiles, we display the difference in magnitude, 
r e , (lM) e , and the FP parameter (logr e — /31og(J) e , see text) between input and recovered values. Points enclosed by triangles have neighbors 
within 1 r e of the primary galaxy being fit. 
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Fig. 6. — The open histogram displays the range of measured Sersic parameter n values found for our real galaxies in clusters RX J0152.7- 
1357, RX J1226. 9+3332, and RX J0142. 0+2131. The solid histogram shows the Sersic n values for only those galaxies which we use to study 
the cluster FP. 




Fig. 7. — Differences between input parameters and values recovered from r 1 / 4 law profile fitting, but for only those galaxies created with n 
> 2 Sersic profiles, within the range of the real galaxies used in the FP. Symbols are the same as in Figure[5] The results for the high redshift 
sample are on the left, lower redshift sample on the right. 
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Fig. 10. — Differences in r 1 / 4 law profile recovered parameter values for n > 2 galaxies as a function of nearest neighbor distance. Points 
enclosed by triangles have neighbors within 1 r e . The average difference in recovered - input values and standard deviation are provided for 
the ranges d < 2, 2 < d < 6, and d > 6 arcsec. 
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Fig. 11. — Comparison of the derived Fundamental Plane parameter for galaxies observed in multiple images. Plotted are the differences 
in repeat measurements from overlapping pairs of images for both RX J0152. 7-1367 and RX J1226. 9+3332. Average differences for each pair 
are shown by the dashed lines. 




Fig. 12. — Offset in magnitude calibration between this work and Blakcslec ct al. (2006) (black). Galaxies with fitted Sersic n values 
differing by less than 0.1 are circled. Boxed points correspond to galaxies found by us to have best fit n > 4.2. The best straight line fit to 
these points is shown. In gray, differences between our simulated galaxy input and r 1 / 4 law recovered parameters. 
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Fig. 13. — Differences between this work and Blakcslce ct al. (2006) in measured structural parameter values (in all cases (this work - 
Blakeslee)) after correction for magnitude zero point differences. Points corresponding to galaxies best fitted with n > 4.2 are boxed. Galaxies 
with measured Sersic n differing by less than O.f are circled. 




Fig. 14. — Top two panels: differences in measured magnitude and FPP from Sersic and r 1 / 4 law fits as a function of Sersic n for all real 
sample galaxies in RX J0152. 7-1357 and RX J1226. 9+3332. Next two panels: the difference in measured magnitude and surface brightness 
vs. measured r e between Sersic and r 1 / 4 law fits. Objects best fit with Sersic n > 2.0 are denoted by large symbols. We show a best straight 
line fit with slope 3.6 to the points in the bottom panel. 
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Fig, 15. — Sersic index vs SExtractor mea sured total magnitude for galaxies in our RX J0152. 7-1357 and RX J1226. 9+3332 FP samples 
Jj0reensen et alj|2006l , [20071 ; IBarr et al.ll2006l) . 



